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I.  INTRODUCTION 

A.        BACKGROUND 

Since  the  1990s,  radar  imaging  and  multi-dimensional  radar  systems  have  become 

widely  used  in  many  areas,  both  military  and  commercial.  This  came  about  due  to  the 
advent  of  more  advanced  waveform  generation  and  signal  processing  techniques.  The 
requirement  for  high  resolution  radars  in  military  applications  has  always  been  a  top 
priority  in  radar  system  design.  However,  there  is  a  fundamental  tradeoff  between 
resolution  and  the  radar  system's  bandwidth.  The  wide  bandwidth  required  for  high 
resolution  is  difficult  to  achieve.  The  stepped  frequency  waveform  is  an  alternative  to 
building  wide  bandwidth  devices. 

Essentially,  the  stepped  frequency  radar  employs  inter-pulse  pulse  compression, 
in  which  the  carrier  frequency  of  successive  pulses  is  changed  by  a  fixed  frequency  step, 
Af.  Due  to  this  technique,  a  wide  bandwidth  is  inherently  achieved  over  several  pulses 
instead  of  within  a  single  pulse.  However,  its  instantaneous  bandwidth  remains  small, 
and  this  feature  is  important  because  of  the  simpler  hardware  involved  in  the  receiver  for 
a  smaller  instantaneous  bandwidth.  The  effective  stepped  frequency  waveform  (SFWF) 
bandwidth  is  the  product  of  the  number  of  coherently  integrated  pulses,  N,  and  the 
frequency  step  size,  Af. 

The  returns  of  the  N  pulses  at  different  carrier  frequencies  form  the  frequency 
spectrum  of  the  target  reflectivity  and  are  coherently  integrated  by  the  inverse  fast 
Fourier  Transform  (IFFT)  to  give  the  high  resolution  range  (HRR)  profile.    The  HRR 


profile  does  not  give  the  absolute  range  position  of  the  target  but  gives  the  relative 
position  of  each  scatterer  within  the  conventional  range  bin. 

The  key  advantage  of  the  stepped  frequency  implementation  as  compared  to  other 
pulse  compression  techniques  is  its  high  range  resolution  (on  the  order  of  a  fraction  of  a 
meter  is  achievable)  without  compromising  its  instantaneous  bandwidth,  or  target 
detection  range.  It  is  easily  implemented  by  adding  a  frequency  synthesizer  to  a 
conventional  narrowband  receiver  and  performing  the  necessary  signal  processing  on  the 
received  signal. 

The  two  main  limitations  of  the  stepped  frequency  radar  implementation  are  as 
follows: 

a.  The  stepped  frequency  radar  cannot  directly  measure  the  velocity  of 
moving  targets  based  on  the  doppler  shift.  This  is  due  to  the  frequency 
shift  in  the  carrier  during  transmission. 

b.  The  stepped  frequency  implementation  requires  additional  time  to  transmit 
and  receive  N  pulses. 

Inverse  Synthetic  Aperture  Radar  (ISAR)  is  used  to  image  targets  in  the  cross- 
range  dimension  due  to  their  rotational  motion.  The  ISAR,  like  the  SFWF,  works  on  the 
principle  of  doppler  shifts  due  to  target  rotation,  which  causes  a  scatterer  on  the  target  to 
change  its  relative  position  from  the  radar  over  time.  It  can  be  seen  that  a  scatterer  located 
in  the  center  of  the  rotating  platform  will  therefore  experience  no  doppler  shifts,  hence 
registering  zero  frequency  shift  over  time.  However,  this  processing  requires  the  phase 
of  the  transmitted  pulse  be  almost  constant  over  the  processing  interval  for  the  Range 
Rate  Interval  (RRI),  which  is  the  duration  of  the  stepped  frequency  waveform. 


Both  of  the  above  techniques  are  also  termed  as  multiple  pulse  processes.  The 
ISAR  technique  requires  a  fixed  frequency  measurement  of  a  non-stationary  target  and 
the  SFWF  requires  a  sweep  frequency  measurement  of  a  stationary  target.  Proper 
selection  of  the  radar  parameters  can  yield  a  radar  design  for  optimum  range  resolutions 
in  both  down-range  and  cross-range  for  a  particular  target  of  interest,  especially  so  for 
radars  with  adaptive  operating  parameters. 

B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  use  two  proven  high  resolution  processing 
techniques,  namely  stepped  frequency  and  ISAR  processing,  to  image  the  simulated 
signal  returns  from  multiple  scatterers  on  an  object.  Applications  include  imaging  for 
target  identification  and  radar  cross  section  (RCS)  measurements.  RCS  measurement  and 
reduction  are  the  primary  applications  of  interest  in  this  research.  Unlike  the  general 
imaging  case  in  which  the  target  and  radar  geometrical  relationships  are  often  non- 
optimum,  the  measurement  application  allows  the  target  motion  to  be  controlled 
precisely.  The  resulting  images  can  isolate  individual  scattering  sources  on  a  target. 
Knowledge  of  the  individual  scatterer's  characteristics  permits  the  development  of  a  RCS 
reduction  strategy. 

C.  RELATED  WORK 

The  theory  of  stepped  frequency  has  been  explored  in  two  previous  theses  [1,2]. 
In  [1],  various  radial  velocity  compensation  techniques  were  explored  to  eliminate  the 
image  spreading  and  shifting  effects,  while  in  [2],  the  effects  of  changing  the  various 


radar  parameters  on  the  stepped  frequency  performance  were  studied.  This  thesis  goes 
one  step  further,  in  that,  the  stepped  frequency  processing  is  combined  with  the  ISAR 
processing  to  form  images  from  the  signal  returns  of  targets  [3]. 

D.        ORGANIZATION 

This  thesis  starts  out  with  an  introduction  of  the  two  high  resolution  processing 
techniques.  Chapter  II  reviews  the  stepped  frequency  principles  [4]  and  Chapter  III 
reviews  the  ISAR  processing  technique  [4,5].  The  considerations  for  selection  of  the 
various  radar  parameters  to  obtain  the  required  radar  performances  in  detecting  and 
imaging  an  aircraft  target  are  also  covered  in  these  two  chapters.  Chapter  IV  combines 
these  two  techniques  in  the  generation  of  a  range-doppler  profile  using  the  Fast  Fourier 
Transform  (FFT)  and  discusses  the  effects  of  the  target's  radial  and  angular  motion  on  the 
processed  target  image,  which  results  in  an  optimal  tracking  window  for  the  selected 
radar  parameters.  This  is  followed  by  the  design  of  the  radar  parameters  to  detect  a  15  m 
aircraft  target  and  the  simulation  of  the  aircraft  image  (both  a  simple  and  complex 
aircraft)  in  Chapters  V  and  VI  respectively.  Chapter  VI  also  shows  the  simulated  effects 
of  range  walk  and  range  spread  on  the  image  of  each  scatterer  when  the  radar  conditions 
for  the  optimal  tracking  window  are  not  met.  Chapter  VII  presents  the  conclusions  an 
suggestions  for  future  work. 


II.        STEPPED  FREQUENCY  PRINCIPLES 

A.         Block  Diagram 

The  hardware  implementation  of  the  stepped  frequency  radar  is  shown  in  Figure 
2.1.  The  coherent  stepped  frequency  synthesizer  generates  an  output  which  is  stepped  in 
frequency  from  pulse  to  pulse  by  a  fixed  frequency  step  size. 
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Figure  2.1  Block  Diagram  of  Step  Frequency  Radar  Receiver. 


On  the  transmitter  side,  the  output  of  the  coherent  oscillator  (COHO)  is  mixed 
with  the  output  of  the  stepped  frequency  synthesizer  (up-conversion).  The  sum  of  these 
two  frequencies  is  up-converted  to  the  final  transmitted  frequency  after  mixing  with  the 
stable  local  oscillator  (STALO).  The  frequency  is  then  pulse-modulated  and  power 
amplified  before  transmission  into  space. 

The  frequency  of  each  transmitted  pulse  is  comprised  of  three  components:  the 
fixed  intermediate  frequency  of  the  coherent  oscillator,  fixed  radio  frequency  of  the  local 
oscillator  and  the  variable  frequency  of  the  stepped  frequency  synthesizer. 

On  the  receiver  side,  the  return  radio-frequency  target  signal  is  amplified  and 
down-converted  by  first  mixing  it  with  the  output  of  the  stable  local  oscillator.  The 
resulting  signal  is  further  down-converted  to  the  intermediate  frequency  (IF)  by  mixing  it 
with  the  output  of  the  stepped  frequency  synthesizer. 

The  stepped  frequency  synthesizer  is  synchronized  with  the  transmitter  to  allow 
the  transmitted  carrier  frequency  and  the  output  of  the  stepped  frequency  synthesizer  to 
use  the  same  stepped  frequency  in  the  same  pulse  repetition  interval.  This  permits  the 
multiple-time-around  clutter  to  be  filtered  out,  as  it  will  lie  outside  the  passband  of  the 
intermediate  frequency  amplifier. 

The  output  signal  obtained  after  the  second  mixing  is  the  intermediate-frequency 
signal  which  is  fed  into  the  coherent  phase  detector  through  an  intermediate-frequency 
amplifier  with  a  bandwidth  approximately  equal  to  the  inverse  of  the  pulse  width, 
centered  on  the  frequency  of  the  coherent  oscillator. 


B.         Waveform 

The  stepped  frequency  waveform  is  the  interpulse  version  of  the  pulse 
compression  technique  employed  in  conventional  high  resolution  radars.  A  sequence  of  N 
pulses  at  a  fixed  pulse  repetition  frequency  (PRF)  but  different  carrier  frequency  is 
transmitted  by  the  stepped  frequency  radar.  Each  set  of  N  pulses  is  called  a  coherent 
processing  interval.  Figure  2.2  shows  a  stepped  frequency  waveform,  and  Figure  2.3 
shows  its  corresponding  frequency  spectrum.  Note  that  the  frequency  of  each  pulse  in  the 
sequence  is  linearly  increased  from  pulse  to  pulse  by  a  fixed  increment,  with  the 
frequency  of  the  first  pulse  being  the  nominal  carrier  frequency,  fo,  and  the  subsequent 
pulses  having  a  difference  of  Af  (frequency  step  size)  from  the  preceding  pulse.  Thus, 
for  the  n-th  transmitted  pulse  in  the  sequence,  with  n=0,l,2,  ....N-l,  the  carrier  frequency 
is  given  by 

f„  =  fo  +  nAf.  (2.1) 

If  the  transmitted  pulsed  signal  for  the  n-th  pulse  is 

si(t)  =  Aicos[27i(f0  +  nAf)t],  (2.2) 

then  the  target  signal  received  after  a  time  delay  of  2R/c  for  a  target  range  of  R  is  given 
by 

s2(t)  =  A2cos[27t(f0  +  nAf)(t-2R/c)],  (2.3) 

where  c  is  the  speed  of  light. 
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Figure  2.2  Pulse  Sequence  of  Step  Frequency  Waveform. 
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Figure  2.3  Frequency  Spectrum  of  Step  Frequency  Waveform. 


The  output  from  the  coherent  phase  detector  is  modeled  as  the  product  of  the 
received  signal  with  the  reference  signal  followed  by  a  lowpass  filter.  The  output  of  the 
mixer  of  the  coherent  phase  detector  is 


2R 

si(t)s2(t)  =  AiA2cos[  27i(f0  +  nAf)t]  cos[27i(f0  +  nAf)(t )  ] 

c 

=  ^^  {cos[  27i(fo  +  nAf)t  -  27i(f0  +  nAf)(t-— )  ]  +  cos[  27i(f0  +  nAf)  —  ]}. 
2  c  c 

(2.4) 

The  first  component  of  Equation  2.4,  which  is  the  high  frequency  term,  is  filtered  out  by 

the  low-pass  filter.  Therefore,  the  output  of  the  coherent  phase  detector  is 


A1A2       _  _    ,.         A_  2R  _ 

s(t)=  cos[  27i(f0  +  nAf) ] 

2  c 


=  Acos[  27i(f0  +  nAf)  —  ].  (2.5) 

c 


The  output  of  the  in-phase  and  quadrature  (I-Q)  detector  can  be  represented  as 

s(t)  =  Ae-J(t,"  =  I1+jQi,  (2.6) 


2R 

where  <{>n  =  27r(f0  +  nAf) 

c 


If  the  target  moves  with  a  constant  radial  velocity,  v,  the  range  changes  with  each 


pulse  as 


Rn  =  R0  +  vnT, 


(2.7) 


where  T  is  the  pulse  repetition  interval  (PRI)  of  the  sequence  of  pulses  in  seconds.     From 
Equations  2.6  and  2.7,  the  phase  can  be  simplified  to 


<\>n  =  27i(f0  +  nAf)2(R0  +  vnT)/c, 


(2.8) 


or  when  expanded  out, 


<}>n  = 


4Ttf0R0      47iAfR0(n-l)T 
+ 
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Tc 
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2r^2 


47tvf0(n-l)T     47iAfv(n-irT 
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Tc 

li 

4*  term 


(2.9) 


The  first  term  represents  a  constant  phase  shift  and  does  not  figure  into  the  stepped 
frequency  implementation.  The  second  term  represents  the  frequency  shift  due  to  the 
frequency  rate  of  change  multiplied  by  the  round-trip  delay  time.  This  is  also  referred  to 
as  the  induced  phase  shift.  This  term  gives  the  desired  range  resolution  of  the  stepped 
frequency  waveform.  The  third  term  of  Equation  2.9  gives  the  doppler  frequency  due  to 
the  radial  motion  of  the  target.  This  component,  nevertheless,  will  add  on  to  the  phase 
component  of  the  2nd  term.  It  will  be  interpreted  as  though  there  is  a  range  movement  of 
the  target.  The  fourth  term  in  Equation  2.9  is  the  frequency  shift  from  pulse  to  pulse  and 
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occurs  due  to  the  relative  velocity  between  target  and  radar.  Like  the  third  term,  it  will 
result  in  target  spreading  in  range,  attenuation  of  target  amplitude  due  to  the  spreading, 
and  shifting  (range  walking)  of  the  target  signal  in  the  range  profile.  Hence,  contributions 
from  the  third  and  fourth  terms  will  corrupt  the  HRR  profile. 

C.  Key  Parameters 

The  following  are  the  key  parameters  used  in  the  stepped  frequency  radar: 

1.         Instantaneous  and  Effective  Bandwidth 

The  instantaneous  bandwidth  of  the  waveform  is  the  same  as  that  of  a  single  pulse, 
and  is  given  by  the  inverse  of  the  pulse  width,  t,  as 

Binst=  i.  (2.10) 

T 

The  effective  bandwidth  is  correspondingly  given  as 

Beff=NAf,  (2.11) 

and  the  design  of  the  stepped  frequency  waveform  is  such  that 

Beff>Bjnst.  (2.12) 
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2.         Downrange  Resolution 

The  downrange  (range  of  target  from  radar)  resolution  of  a  conventional  radar  is 


ex 

^^conventional    =  ~Z~  ■  (2. 13) 


With  the  stepped  frequency  implementation,  the  downrange  resolution  can  be  made  much 
smaller,  as  follows: 


c  c 

Ard  =  =  .  (2.14) 

2Beff  2NAf 


From  Equation  2.14,  the  downrange  resolution  of  the  stepped  frequency  waveform  can  be 
increased  by  either  increasing  the  number  of  pulses  in  the  sequence  or  the  frequency  step 
size.  Alternatively,  the  increase  in  the  downrange  resolution  can  be  seen  as  the  original 
range  resolution  cell  being  broken  up  into  N  sub-cells.  The  FFT  output  of  the  N  samples 
will  generate  a  high  resolution  range  profile,  with  resolution 


Ru           c                 ct/2 
Ard  =—   =    =    ,  (2.15) 

N  2NAf  NxAf 


where  Ru  is  the  unambiguous  range  window. 
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A  simplified  pictorial  representation  of  Ard,COnventionai  ,  Ard,  and  Ru  is  given 
in  Figure  2.4.  From  Equation  2.15,  it  can  be  seen  that  the  product  of  pulse  width  and 
frequency  step  size  determines  the  downrange  resolution  of  the  stepped  frequency  radar. 


Figure  2.4  Range  Resolution  with  Step  Frequency  Radar. 

3.  Downrange  Ambiguities 

Downrange  ambiguities  exist  for  the  stepped  frequency  radar.  The  maximum 
unambiguous  downrange  is  associated  with  the  SFWF  process  and  can  be  determined 
from  the  phase  (J)n  of  the  received  echo  for  frequency  step  n  from  a  scatterer  at  range  R0, 
given  by  the  2nd  term  in  the  target  signal  phase  Equation  2.9 


_      .     Af2Ro.     ,._ 

<&  =  2n (n-l)T. 

Tc 


(2.16) 
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The  corresponding  frequency  shift  can  be  determined  from  the  rate  of  phase 


change  as 


Af  2R0 

f=   .  (2.17) 

Tc 


By  equating  the  frequency  shift  to  the  PRF,  fp 


fp  =  l/T=^.  (2.18) 

Tc 


Simplifying,  the  downrange  ambiguity  window,  Ru,  which  is  equivalent  to  the  maximum 
downrange  image  dimension,  Wd,  gives 


Ru=—  •  (2.19) 

2Af 


For  a  frequency  step  size,  Af,  of  10MHz,  the  unambiguous  range  window  is  15m. 
Scatterers  located  at  10m,  25m,  625m  will  all  appear  to  be  at  the  same  range.  The 
stepped  frequency  processing  cannot,  therefore,  provide  absolute  range  to  scatterer,  but 
rather  relative  range  position  within  the  received  pulse.  The  normal  time  delay 
measurements  (range  gating)  will  have  to  be  used  to  compute  the  actual  range  to  the 
target  under  range  track.  Therefore,  the  range  of  each  scatterer  is  computed  from  the 
combination  of  information  from  range  gating  and  relative  range  from  stepped  frequency 
processing.  The 
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absolute  range  to  a  scatterer  is  given  as 


R=  R0  +  n 


2NAf 


(2.20) 


where  Ro  is  the  absolute  range  of  the  scatterer  from  the  beginning  of  the  unambiguous 
range  window.  The  start  range  of  the  unambiguous  range  window  is  computed  as  a 
integer  (n)  multiple  of  c/2Af.  For  example,  with  n  =  60,  the  range  gate  opens  from  900m 
to  915m.  Within  that  interval,  the  smallest  resolution  is  given  in  Equation  2.15. 

4.  Product    xAf 

The  high  range  resolution  (HRR)  profile  of  the  stepped  frequency  radar  depends 
to  a  large  extent  on  the  product  of  the  pulse  width  and  frequency  step  size,  xAf.  There  are 
three  cases: 

a)  xAf  =  1.  The  unambiguous  downrange  window,  Ru,  equals  the 
original  range  bin  size  of  cx/2.  This  represents  the  minimum 
required  sampling  rate  in  the  frequency  domain  for  a  given  pulse 
width.  The  disadvantage  of  this  condition  is  that  the  wraparound 
of  scatterer  peaks  caused  by  slight  changes  in  the  target  range  is 
difficult  to  notice. 

b)  xAf  <  1.  The  original  range  bin  occupies  only  part  of  the 
unambiguous  downrange  window.  This  implies  oversampling  in 
the  frequency  domain  and  will  result  in  a  loss  in  downrange 
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resolution  in  the  range  profile.    The  advantage  of  this  is  that  the 
unused  part  of  the  downrange  profile  serves  as  a  clutter-free  zone 


to  detect  moving  targets. 


c)  xAf  >  1 .  The  unambiguous  downrange  window  is  smaller  than  the 

original  range  bin.  This  represents  undersampling  in  the  frequency 
domain.  Folding  over  of  targets  will  result  and  therefore  this 
condition  should  be  avoided. 

5.  Selection  of  Frequency  Step  Size 

The  frequency  step  size  is  important  in  determining  the  size  of  the  ambiguous 
downrange  window  and  the  range  resolution  under  stepped  frequency  implementation. 
Generally,  the  frequency  step  size  is  chosen  to  be  approximately  half  the  inverse  of  the 
pulse-width  [2].  In  cases  where  there  are  no  background  scatterers  outside  the  target,  like 
working  in  an  anechoic  chamber,  the  target  maximum  extent,  E,  plays  a  part  in 
influencing  the  frequency  step  size,  as  follows: 


Af<-^.  (2.21) 

2E 


From  the  above  relationship,  the  frequency  step  size  varies  inversely  as  the  target 
extent,  E.   In  addition,  the  pulse  width  needs  to  be  large  enough  to  encompass  the  entire 
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target 


2E 

x> — .  (2.22) 

c 


Combining  Equations  2.18  and  2.19,  we  obtain  the  requirement  that 


1 
Af<— .  (2.23) 

T 
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III.       INVERSE  SYNTHETIC  ARRAY  RADAR  (ISAR)  PROCESSING 

A.  ISAR  THEORY 

ISAR  provides  a  processing  technique  capable  of  high  crossrange  resolution  due 
to  target  rotation.  This  method  assumes  that  velocity  compensation  has  been  done  to 
eliminate  target-to-radar  radial  velocity  component,  which,  as  explained  in  the  previous 
chapter,  has  the  effect  of  range  bin  shift  (or  range  walk)  as  well  as  target  spreading  within 
a  bin.  We  recall  that  the  SFWF  generates  a  frequency  shift  from  a  known  frequency 
change  on  a  pulse-to-pulse  basis.  It  requires  that  the  target  be  stationary  or  velocity- 
compensated  during  the  RRI.  The  ISAR  processing  measures  the  dominant  scatterers' 
locations  on  a  target  arising  from  a  shift  in  doppler  frequency  due  to  changing  scatterers' 
radial  velocities  as  the  target  rotates. 

B.  ISAR  CROSSRANGE  CALCULATION 

The  crossrange  resolution  obtained  from  the  ISAR  processing  is  dependent 
primarily  on  the  doppler  shift  of  the  various  scatterers  located  on  the  target.  This  is 
illustrated  by  Figure  3.1.  For  this  example,  there  is  no  translational  motion  between  the 
radar  and  the  target  and  therefore  the  depression  angle  is  zero;  that  is,  all  scatterers  on  the 
target  and  the  radar  are  in  the  same  horizontal  plane.  Also,  the  target  rotates  at  a  constant 
angular  rotation  speed  of  u)  rad/s  about  the  axis  that  is  perpendicular  to  the  horizontal 
plane.  From  the  figure,  it  can  be  seen  that  a  point  scatterer  at  a  crossrange  of  rc  from 
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the  center  of  rotation  will  have  a  instantaneous  tangential  velocity,  v,  given  by 


v=cor. 


(3.1) 


O  Scatterer  on  a  target 
•  Center  of  rotation 


Radar 


V2    -4- 


LOS 


($t=AQ. 


Figure  3.1  ISAR  Crossrange  Resolution  of  Two  Scatterers. 


The  corresponding  doppler  frequency  can  be  computed  from 


fd  =  2corc  - , 
c 


(3.2) 


where  fc  is  the  carrier  frequency  and  c  the  speed  of  light  at  3  x  10    m/s.  For  the  case  of  2 
point  scatterers  located  in  the  same  downrange  cell,  but  separated  in  crossrange  by  Arc, 
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the  tangential  velocity  difference  is  a)Arc,  and  the  corresponding  frequency  difference  can 
be  computed  as 


Afc  =  2co  ^.  (3.3) 

c 


The   doppler  resolution    is   approximately   equal   to   the   reciprocal    of  the   coherent 
integration  time, 


Afd  =  1/Tj  = ! ,  (3.4) 

MxPRI 


where  M  is  the  total  number  of  pulses  processed  for  doppler  resolution.    For  a  doppler 
resolution  of  Afd,  the  corresponding  crossrange  resolution  is 


cAfd 
Arc= -,  (3.5) 

2cofc 


and  from  Equations  3.3,  3.4  and  3.5,  we  obtain  another  expression  for  the  crossrange 
resolution: 


Arc=  -^-     =    —  =    — ,  (3.6) 

2(i)Tifc  2a)Ti         2\j/ 
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where  X  is  the  wavelength  and  \\f  the  total  angle  moved  during  the  integration  time.  As 
was  done  for  the  SFWF,  the  Fast  Fourier  transform  (FFT)  is  used  to  convert  the  M  time 
samples  into  a  frequency  spectrum. 

C.         MAXIMUM  UNFOCUSSED  INTEGRATION  ANGLE 

There  is  the  maximum  angle  for  the  target  rotation,  depending  on  the  target 
maximum  length,  whereby  phase  deviation  occurs  which  causes  defocusing  of  the  target 
image  in  the  crossrange  after  processing.  The  defocusing  occurs  because  of  the  change  in 
distance  between  the  radar  and  the  scatterer  as  the  target  rotates.  This  distance  can  be 
converted  to  an  equivalent  phase  difference.  The  rule  of  thumb  for  the  maximum 
allowable  two-way  phase  deviation  before  defocusing  occurs  is  given  as  xt/8  radians  [2]. 
The  maximum  integration  angle  before  defocusing  occurs  is 


¥max=-j-.  0.7) 

2  Vr 


E 

The  maximum  target  radius,  r  (=  —  if  the  center  of  the  target  is  the  center  of  rotation), 

2 


before  defocusing  occurs,  is  expressed  as 


2(Arc)2 
r<    .  (3.8) 
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If  the  criterion  for  \j/max  is  violated,  then  focusing  is  required.   Essentially,  this  is 
done  by  subtracting  or  adding  the  two-way  path  difference,  hence  phase  difference,  as  the 


scatterer  rotates  around  the  center  of  the  target. 


D.        CROSSRANGE  SAMPLING  CRITERIA 

The  selection  of  PRF  is  determined  by  the  Nyquist  sampling  requirement  to 
prevent  aliasing.  Nyquist  sampling  requires  that  the  signal  be  sampled  at  twice  the 
highest  expected  doppler  frequency  component.  The  highest  doppler  frequency  is 
obtained  from  the  scatterer  that  is  furthest  away  from  the  center  of  rotation,  and  is  given 
by 


4fc 

r  Rreffective  =  fp.effective  —  •^ld(max)  =  tOrC(max),  {->•?) 

C 


where  rc(max)  is  the  distance  of  the  furthest  scatterer  from  the  center  of  rotation.  The 
effective  PRF  for  a  stepped  frequency  radar  can  also  given  be  given  as 

PRFeffective=PRF/N,  (3.10) 

where  N  is  the  number  of  frequency  steps  used.  This  is  the  unambiguous  doppler 
window.  Therefore,  the  maximum  and  minimum  scatterer  doppler  frequencies  that  will 
not  cause  aliasing  are  from  +PRFeffective/2  to  -PRFeffective/2  Hz.  The  corresponding 
unambiguous  crossrange  window  is  given  as 
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C  P Kreffecuve  MA, 

Wc  =  2rmax=  =    ,  (3.11) 

2cofc  2coTi 


where  Tj    is  the  total  coherent  integration  time  for  the  M  ISAR  waveforms  (fixed 
frequency  time  samples). 
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IV.       RANGE-DOPPLER  PROCESSING 

A.         UNAMBIGUOUS  WINDOW  AND  RESOLUTION 

The  two-dimensional  image  of  a  target  can  be  generated  from  the  stepped 
frequency  and  ISAR  processing.  The  two  image  axes  are  downrange  and  doppler 
(crossrange).  The  downrange  and  crossrange  resolutions  are  given  in  Equations  2.15  and 
3.6,  while  the  downrange  and  crossrange  windows  are  given  in  Equations  2.21  and  3.10, 
respectively.  Both  the  downrange  and  doppler  (crossrange)  windows  will  combine  to 
give  the  desired  two-dimensional  image  of  a  target. 

In  order  to  obtain  these  two  dimensions,  two  transforms  are  required.  For  the 
SFWF,  N  pulses  are  transmitted  in  one  waveform  and  during  these  N  pulses,  it  is 
assumed  that  there  is  only  a  very  small  viewing  angle  change.  For  each  of  the  N  pulses,  a 
range  profile  can  be  generated.  Altogether,  there  are  M  such  range  profiles,  which  can  be 
transformed  by  a  discrete  Fourier  transform  (DFT)  to  obtain  the  target  doppler 
information  for  each  range  cell.  The  N  pulses  and  M  profiles  will  therefore  define  the 
downrange  and  crossrange  range  windows. 

Figure  4.1  illustrates  the  SFWF-ISAR  waveform  relationship.  There  are  M  (=4) 
stepped  frequency  waveforms,  each  consisting  of  N  (=16)  stepped  frequency  pulses.  If 
the  RRI  (=  N  x  PRI)  is  equal  to  the  dwell  time,  ie,  there  is  no  idle  time,  then  the  total 
integration  time  T;  is 

Ti  =  M  x  RRI 

=  MxNxPRI.  (4.1) 
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Figure  4.1  Frequencies  of  Successive  Pulses  in  the  SFWF.  Each  SFWF  is 
Processed  to  Form  One  of  M  HRR  Profiles  Used  In  Generating  The  ISAR  Profiles. 


B.         FAST  FOURIER  TRANSFORM  (FFT) 

To  illustrate  how  the  FFT  is  used  to  generate  the  two-dimensional  image  of  the 
target,  the  returned  M  x  N  data  pulses  collected  are  organised  into  a  2-D  array  as  shown 
in  Figure  4.2.  The  received  frequency  signatures  of  the  M  bursts  are  the  time  history 
series  of  the  target's  reflectivity  at  each  discrete  frequency.  The  radar  processor  then 
uses  the  frequency  signatures  as  the  raw  data  to  perform  range  and  doppler  processing. 
For  the  stepped  frequency  signals,  the  range  processing  performs  an  N-point  1-D  inverse 
Fourier  transform  (1-EDFT)  for  each  of  the  M  received  frequency  signatures.  Therefore, 
M  range  profiles  (ie,  the  distribution  of  the  target  reflectivities  in  range),  each  containing 
N  range  cells,  can  be  obtained.  For  each  range  cell,  the  M  range  profiles  constitute  a  new 
time  series,  which  is  sampled  at  baseband  to  form  M  in-phase  (I  channel)  and  M 
quadrature-phase  (Q-channel)  data. 
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Figure  4.2  Illustration  of  a  Step  Frequency  ISAR  Imaging  of  a  Rotating  Target  (after  [3]). 
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The  1-D  FFT  (1-DFFT)  is  then  used  for  the  time  history  series  to  generate  an  M-point 
doppler  spectrum,  or  doppler  profile.  Hence,  a  N  x  M  image  is  formed  by  combining  the 
N  doppler  spectra  at  N  range  cells.  Figure  4.2  illustrates  the  various  steps  in  the 
processing. 

C.         TARGET  RADIAL  (TRANSLATIONAL)  VELOCITY  EFFECTS 

In  this  thesis,  it  is  assumed  that  only  rotational  motion  of  the  aircraft  target  exists; 
no  translational  velocity  is  modeled  or  if  it  is  present,  is  assumed  to  be  compensated. 
However,  for  completeness,  it  is  discussed  here. 

Translational  velocity  causes  the  effects  described  in  the  following  two  sections. 

1.  Range  Offset 

Also  known  as  range  spreading,  this  causes  the  range  of  each  scatterer  to  spread 
from  its  actual  location  to  adjacent  range  cells,  thus  providing  a  distorted  and  undesirable 
range  profile.  This  spreading  needs  to  be  compensated.  The  formula  that  relates  the 
range  offset  or  number  of  bins  spread,  P,  to  the  radial  (translational)  velocity  is  given  as 
follows: 


vrNT 
P=    —  (4.2) 

Ar 


where  NT  =  total  duration  of  the  stepped  frequency  waveform,  Ar  is  the  effective  range 
resolution,  and  vr  the  radial  velocity  of  target. 
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2.  Range  Walk 

This  has  the  effect  of  shifting  the  range  profile  by  L  bins,  resulting  in  an 
inaccurate  position  processing  of  individual  scatterers.  The  formula  that  relates  this  range 
bin  shift  to  the  radial  velocity  is 


L=i£p=     Si*L,  (4.3) 


Beff  f, 


P 


where  fd  is  the  doppler  due  to  uncompensated  radial  velocity,  ie, 


fc'Vr 

fd  =  ,  (4.4) 

BeffAr 


and  fc'  is  the  nominal  transmitter  frequency  (center  frequency) 

fc'  =       fo+Beff.  (4.5) 

D.         TARGET  ANGULAR  VELOCITY  EFFECTS 

There  are  two  effects  due  to  the  rotational  motion  of  the  target:  (1)  range-doppler 
interaction  and  (2)  cell  migration. 

1.  Range-Doppler  Interaction 

This  will  cause  range  smearing  to  the  range  profile.  As  a  rule  of  thumb  [1],    in 

order  to  avoid  range  smearing,  the  rotation  angle  of  the  target,  G),  should  satisfy 
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1  cfp 

w<v .  (4.6) 

2  ENfc 


The  factor  of  Vi  is  chosen  as  a  rule  of  thumb. 

2.  Cell  Migration 

Distortion  of  the  range  profile  occurs  due  to  the  quadratic  phase  variation  caused 
by  the  circular  motion  of  the  scatterers  as  they  rotate.  Another  way  of  looking  at  this  is 
that  the  radial  velocity  of  the  scatterer  is  not  constant  throughout  the  integration  time,  Tj, 
during  the  rotation.  When  used  with  the  FFT,  movement  of  the  scatterers  across  range 
cells  occur  and  blurring  results.  The  movement  across  range  cells  can  occur  in  both 
crossrange  and  downrange. 

In  order  to  prevent  cell  migration  from  occurring,  the  selection  of  radar 
parameters  is  important.  The  occurrence  of  image  blurring  can  be  represented  by  an 
indicator,  Nd  (for  downrange)  or  Nc  (for  crossrange)  as  follows: 


E 

_2Arc,d_ 

c 

.fc. 

1 

_2Arc.d_ 

NCd  = —     •  (47) 

T.2Arc,dJLfcjL2Arc.dJ  ^,} 

Nc  and  Nd  represent  the  number  of  range  bin  shifts  the  image  experiences  in  the 
crossrange  and  downrange  respectively.  To  prevent  blurring  of  the  image,  Ncd  should  be 
less  than  unity.  A  factor  of  V2  is  used  as  a  rule  of  thumb.  From  Equations  2.14,  2.19,  3.6, 
and  3.11,  the  values  of  Nc  and  Nd  determine  the  optimal  crossrange  and  downrange 
window.  A  scatterer  located  within  this  window  will  experience  no  range  bin  shifts. 
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V.         SIMULATION  DESIGN 

A.         APPROACH 

In  this  thesis,  various  point  scatterer  models  of  aircraft  will  be  used  as  the  targets 
from  which  signals  will  be  received  by  the  radar  and  processed  using  stepped  frequency 
and  ISAR  techniques.  Matlab  will  be  used  to  simulate  the  target  returns  and  provide  the 
processing  tools  for  the  BFFT  and  DFFT.  The  approach  adopted  in  the  Matlab 
programming  is  as  follows: 

1.  The  aircraft  scatterers  are  computed  in  cartesian  co-ordinates  manually  on 
paper  or  generated  using  a  Computer-Aided  program.  The  radar  is  located 
at  the  origin. 

2.  A  rotation  speed  at  co  radians  per  second  is  selected  so  that  it  meets  the 
Nyquist  sampling  criteria. 

3.  The  various  radar  parameters  are  selected  and  verified  to  meet  both 
stepped  frequency  and  ISAR  processing  requirements. 

4.  The  radial  velocities,  hence  doppler  frequencies,  of  the  scatterers,  as  seen 
by  the  radar,  are  computed. 

5.  Amplitudes  for  the  scatterers  are  varied  according  to  their  orientation  to 
the  radar  to  simulate  a  realistic  target.  For  example,  those  scatterers  on  the 
leading  edge  of  an  aircraft  are  assigned  a  higher  radar  cross  section 
compared  to  the  rest  of  the  scatterers  on  the  trailing  edge  of  the  aircraft. 
This  applies  only  to  the  complex  target. 
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6.  The  radar  transmits  M  pulse  bursts,  each  with  N  pulses,  stepped  in 
frequency.  It  is  assumed  that,  within  each  burst,  the  target's  doppler  is 
constant.  Only  the  doppler  frequency  of  each  of  the  M  bursts  changes. 

7.  The  I  and  Q  components  of  the  target  returns  are  generated  for  each  pulse. 

8.  The  I  and  Q  components  of  the  N  pulses  within  each  burst  are  processed 
using  the  IFFT  to  generate  the  M  HRR  (downrange  information)  profiles. 

9.  The  M  HRR  profiles  obtained  from  step  8  are  further  processed  using  the 
DFFT  to  generate  the  Doppler  profile.  The  doppler  information  is  then 
mapped  into  crossrange  coordinates. 

10.  The  two-dimensional  image  of  the  target  is  mapped  using  the  HRR  and 
doppler  profiles. 

B.         TARGET 

A  simple  wire  model  of  an  aircraft  with  9  scatterers  will  be  modeled  first  in  order 
to  illustrate  the  stepped  frequency  and  ISAR  processing  without  having  to  look  at  the 
enormous  details  arising  from  a  large  number  of  scatterers.  The  dimensions  are  as  shown 
in  Figure  5.1.  The  maximum  target  extent  is  15m  and  it  is  initially  nose-on  to  the  radar. 

In  the  case  of  a  complex  target,  22  scatterers  are  used  to  represent  the  various 
scattering  points  on  the  aircraft.  These  points  realistically  represent  the  dominant 
scatterers  from  an  aircraft.  The  maximum  target  extent  is  10m.  The  dimensions  are  as 
given  in  Figure  5.2.  The  leading  edge  scatterers,  ie,  those  facing  the  radar,  will  be 
assigned  higher  signal  returns.  Note  that,  in  both  cases,  the  target  is  assumed  to  be  in  the 
same  horizontal  plane  as  the  radar  so  that  all  its  z-coordinates  are  zero.    The  rotational 
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speed  of  the  target,  co,  is  assumed  to  be  0.002  rad/s,  which  may  not  typically  represent  the 
yaw  dynamics  of  a  rotating  target.  This  value  of  co  is  chosen  in  order  to  match  the  radar 
parameters  selected  to  avoid  image  distortion.  The  key  point  of  the  ISAR  technique  is 
that  crossranse  information  can  be  extracted  as  long  as  there  is  rotational  motion  of  the 
target.  In  a  real  world,  a  target's  motion  will  have  a  rotating  component  in  which  the 
radar  can  process  to  track  and  therefore  image  the  target's  RCS  in  the  crossrange 
direction. 


25m       27m       30m 


40m 


(downrange) 


Figure  5.1  Model  of  a  Simplified  Aircraft  with  9  Scatterers. 
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Figure  5.2  Model  of  a  Complex  Aircraft  with  22  Scatterers. 
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C.  RADAR  AND  TARGET  PARAMETERS 

The  following  radar  and  target  parameters  were  selected: 

1.  Carrier  frequency,  f0  =  40  GHz 

2.  PRF,  fp  =  10  kHz 

3.  Stepped  frequency  step  size,  Af  =1  MHz 

4.  No.  of  stepped  frequency  steps,  N  =512 

5.  No.  of  burst  for  ISAR  processing,  M  =  100 

6.  Pulse  width,  T  =0.1  jj.s 

7.  Maximum  target  extent,  E  =  15m 
(radar  designed  to  image  small  fighter  aircraft  only) 

8.  Target  rotational  frequency,  u)  =  0.002  rad/s 

9.  Unambiguous  downrange  window,  Ru  (Wd)  =  150  m 

D.  STEPPED       FREQUENCY       PROCESSING       PARAMETERS       AND 
REQUIREMENTS 

The  following  parameters  were  computed  and  verified  to  meet  the  stepped 
frequency  requirements  as  delineated  in  Chapter  II. 

1.  Receiver  instantaneous  bandwidth,  Binst  (=  1/x)  =  10  MHz 

2.  Frequency  step  size,  Af  <  1/x  ?  (106  <  107  ?)  =  YES 

3.  Frequency  step  size,  Af  <  c/(2E)  ?  (106  <  107?)  =  YES 

4.  Effective  signal  bandwidth,  Beff  (=NAf)  =512  MHz 

5.  Conventional  downrange  resolution,  Arconventionai  =  15  m 
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6.  Stepped  frequency  downrange  resolution,  Ar  =  0.293  m 

7.  Conventional  unambiguous  downrange  window,  Ru,COnventionai  =  15km 

8.  Stepped  frequency  unambiguous  downrange  window,  Ru  (=Wd)  =  150  m 

9.  Ru  >E?  (150m  >  15m?)  =  YES 

2F 

10.  Pulsewidth,x>  —  ?    (10"7>10"7?)  =  YES 

c 

E.       ISAR  PROCESSING  PARAMETERS  AND  REQUIREMENTS 

From  the  radar  parameters,  the  ISAR  requirements  were  verified  to  meet  the 
Nyquist  sampling  requirement  to  prevent  aliasing. 

1.  Integration  time,  Ti  =  5.12  s 

2.  doppler  resolution,  Afd  =  0.195  Hz 

3.  Crossrange  resolution,  Arc  =  0.366  m 

4.  Unambiguous  crossrange  window  size,  Wc  =  36.6  m 
(max  target  width  processed) 

5.  Angle  rotated  by  target  during  integration  time,  \j/  =0.01  rad 

6.  Maximum  allowed  rotation  angle,  \j/max  =  0.01 1  rad 

7.  \|/  <i|W?  (0.010  rad  <  0.011  rad?)  =  YES 

8.  Effective  PRF,  fp,effect.ve  =  19.5  Hz 

9.  Twice  doppler  frequency,  2fd(max)  =  17  Hz 

10.  fp.effect.ve  >  2fd(max)  ?  )  (19.5  Hz  >  17  Hz?)  =  YES 
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VI.  SIMULATION  RESULTS  &  DISCUSSION 

The  primary  objective  of  the  simulation  is  to  image  an  aircraft  target  comprised  of 
multiple  scatterers.  As  the  downrange  resolution  of  the  radar  processing  is  dependednt 
on  the  effective  bandwidth,  it  is  important  to  visualize  how  the  effective  bandwidth  of  a 
stepped  frequency  waveform  is  obtained,  and  the  undesirable  effects  of  target  rotation  on 
the  image  quality.  The  first  two  parts  of  this  section  will  focus  on  performing  simulations 
to  illustrate  these  two  concepts.  It  will  be  followed  by  simulating  the  stepped  frequency 
and  ISAR  processing  techniques  to  image  both  simple  and  complex  targets. 

A.         BANDWIDTH  OF  TRANSMITTED  WAVEFORM 

From  Chapter  II,  it  can  be  seen  that  the  high  downrange  resolution  is  obtained 
from  an  wide  effective  bandwidth,  Beff.  Figures  6.1  through  6.3  show  the  simulation 
results  for  the  following  parameters.  The  figures  demonstrate  that  the  effective 
bandwidth  is  NAf. 


a. 

Carrier  frequency,  fc 

=  1GHz 

b. 

Number  of  datapoints  for  FFT 

=  512 

c. 

Stepped  frequency,  Af 

=  5MHz 

d. 

Number  of  pulses,  N 

=  12 

e. 

NAf 

=  60MHz 

f. 

Pulsewidth,  x 

=  0.1us 

o 

©• 

Pulse  repetitive  frequency, 

fP 

=  10kHz 
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Fisure  6. 1    Transmitted  Pulse  Soectrum  after  FFT. 
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Figure  6.2    Transmitted  Pulse  Spectrum  after  FFTSHIFT. 
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Figure  6.3    Expanded  Pulse  Spectrum  after  FFT  and  FFTSHIFT  (Positive  Frequency). 

In  the  simulation,  12  stepped  frequency  pulses  are  transmitted  and  the  signal 
returns  are  processed  using  FFT.  A  characteristic  of  the  function  FFT  is  that  the  positive 
frequencies  lie  on  the  left  side  of  the  spectrum.  FFTSFfTFT  is  used  to  properly  order  the 
spectrum.  The  transmitted  pulses,  FFT  and  FFTSHIFT  processed  spectrums  are  shown  in 
Figure  6.1  and  Figure  6.2.  Figure  6.1  shows  both  the  transmitted  waveform  and  its 
spectrum,  while  Figure  6.2  shows  the  FFTSHIFT  spectra  with  the  x-axis  scaled  as  FFT 
bin  number  and  frequency,  respectively.  By  keeping  only  the  positive  frequency 
spectrum  through  filtering,  the  step  frequency  spectrum  showing  an  increased  bandwidth 
can  be  extracted.  The  simulation  shown  in  Figure  6.3  demonstrates  that  the  stepped 
frequency  waveform  exhibits  an  effective  bandwidth  of  60MHz,  which  is  obtained  from 
the  product  of  the  frequency  step  size,  Af,  and  the  number  of  pulses  transmitted,  N. 
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B. 


TRACKING     WINDOW    AND    EFFECTS     OF    TARGET    ANGULAR 


VELOCITY 

Chapter  IV  discussed  the  effects  of  the  target  angular  velocity  on  both  the  range 
and  doppler  profiles,  namely,   range-doppler  interaction  and  cell  migration.     The  image 
effects  include  shifting  of  range  cells  as  well  as  blurring. 

With  the  radar  and  target  parameters  chosen  for  the  design,  we  can  map  the 
tracking  window  in  which  unambiguous  downrange  and  crossrange  information  can  be 
derived.  The  tracking  window  dimensions  are  determined  by  the  unambiguous 
downrange  (Wd)  and  the  crossrange  window  (Wc),  respectively.  Any  scatterer  that  falls 
outside  this  window  will  wrap  around  and  be  incorrectly  displayed  in  the  tracking 
window.  The  dimensions  of  the  tracking  window  are  shown  in  Figure  6.4. 


Wd(150m)~ 


Wc  =  36.6  m 


Tracking  Window 


Figure  6.4  Tracking  Window. 


The  effect  of  range  wraparound  can  be  shown  from  the  following  case,  whereby 
the  tracking  window  stretches  from  4500  m  to  4650  m  (Ru  =  150m).  A  scatterer  located  at 
4690  m  will  appear  at  40  m  in  the  tracking  window.  This  is  illustrated  in  Figure  6.5. 
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Figure  6.5  HRR  (downrange)  Profile  -  Range  Wraparound. 

If  the  target  rotation  does  not  meet  the  requirements  put  forth  in  Chapter  IV,  the 
result  will  be  a  degradation  of  the  range  and  doppler  (or  crossrange)  profiles,  depending 
on  the  location  of  the  scatterers  from  the  center  of  rotation  and  the  radar-target  line-of- 
sight.  In  the  simulation,  scatterers  were  defined  at  different  locations  from  the  center  of 
rotation.  One  of  the  scatterers  is  always  in  the  center  of  the  rotation  for  comparison 
purposes.  Two  different  scenarios  are  tabulated  in  Table  1.  Note  that  the  threshold  for 
the  downrange  and  crossrange  radial  distances  from  the  center  of  rotation  are  computed 
from  Equation  4.7  to  be  approximately  14  m  and  18  m  respectively.  This  means  that,  for 
this  design,  the  crossrange  window  width  (about  18m  in  radius)  is  the  limiting  factor 
instead  of  the  threshold  distance  from  the  center  of  rotation  to  avoid  spreading  and  range 
walking  effects  of  the  range  profile.  A  scatterer  located  at  the  same  range  to  the  radar  as 
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the  scatterer  located  at  the  center  of  rotation,  but  at  different  crossrange  will  result  in  a 
shifting  of  the  range  profile  in  the  downrange  direction.  This  will  be  illustrated  in  the 
output  graphs  for  scenario  one.  A  similar  explanation  holds  for  scatterers  having  a 
different  downrange  but  with  the  same  crossrange.  The  observation  is  illustrated  in  the 
output  graphs  from  scenario  two. 


Coordinates  (x,y)  of 

Scatterer  1 
(center  of  rotation) 

Coordinates  (x,y) 
of  Scatterer  2 

Coordinates  (x,y) 
of  Scatterer  3 

Coordinates  (x,y) 
of  Scatterer  4 

Scenario 
One 

(4530,0) 

(4540,0) 

(4560,0) 

(4600,0) 

Nd 

- 

<0.5 

<0.5 

<0.5 

Nc 

- 

<0.5 

>0.5 

>0.5 

Scenario 
Two 

(4520,7) 

(4530,5) 

(4530,17) 

Nd 

- 

<0.5 

e0.5 

Nc 

- 

<0.5 

<0.5 

1. 


Table  1  Effects  of  Target  Rotation  -  Two  Scenarios. 


Scenario  One 


In  this  scenario,  the  four  scatterers  are  located  in  a  straight  line  along  the  radar-to- 
target  center  line,  as  shown  in  Figure  6.6a.  Figure  6.6b  shows  the  downrange  profile, 
which  does  not  indicate  any  excessive  (Nc,d)  spreading  effects,  as  expected.  However,  in 
the  crossrange  profile,  scatterers  3  (downrange  separation  of  30m)  and  4  (downrange 
separation  of  70m)  do  experience  undesirable  spreading  and  range  walking  effects.  The 
spreading  and  range  walking  effects  are  shown  in  Figures  6.6c  and  6.6d  respectively. 
This  is  because  both  scatterers  exceeded  the  threshold  radial  distances.  From  these 
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limitations,  the  downrange  tracking  window  has  been  optimized  to  work  in  the  range  of 
about  ±  18m  from  the  target's  center  of  rotation  to  provide  an  undistorted  downrange 
profile. 


downrange 


Scatterer  1 
(center  of— 
rotation) 
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(b) 
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Downrange  in  meters 


Scatterer  on 
enter  of  rotation 


Crossrange  in  meters 


0      -20 
,   <C) 

Figure  6.6  (a)  Location  of  Scatterers,  (b)  Downrange  Profile  (c)  3-D 
View  Depicting  Target  Spreading  and  Range  Walking  Effects. 
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Figure  6.6  (continued)  (d)  Plan  View  Depicting  Target  Spreading  and 


Range 


Walking  Effects. 
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Scenario  Two 


For  this  scenario,  the  crossrange  tracking  window  of  approximately  ±18m  has 
been  noted  in  earlier  paragraphs  as  the  limiting  width  for  an  undistorted  downrange 
profile.  The  Matlab  plots  given  in  Figures  6.7a  through  d  verify  this. 
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Figure  6.7  (a)  Location  of  Scatterers,  (b)  Downrange  Profile, 
(c)  Crossrange  Profile. 
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Figure  6.7  (continued)  (d)  Plan  View  Depicting  Minimal  Target 
Spreading  and  Range  Walking  Effects. 


In  summary,  the  size  of  the  optimized  tracking  window  depends  largely  on  the 
rotational  speed,  the  integration  time  and  the  distance  of  the  scatterer  from  the  center  of 
rotation.  In  the  design,  the  optimized  tracking  window,  whereby  the  target  image 
remains  undistorted  in  both  downrange  and  crossrange  dimensions  due  to  its  rotation,  is 
given  in  Figure  6.8. 


Wc  =  36.6  m 


Wd(150m)~ 


Original  Tracking  Window 


28m 


Optimized 

indow 


Figure  6.8  Optimized  Tracking  Window. 
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C.         IMAGE  OF  SIMPLE  TARGET 

The  image  of  the  simple  target  in  Figure  5.1  is  shown  in  Figures  6.9  and  6.10. 
The  RCS  returns  from  the  9  scatterers  are  clearly  shown  in  the  mesh  plot  of  Figure  6.9. 
All  scatterers  were  defined  within  the  optimized  tracking  window  to  avoid  degradation  of 
the  downrange  and  crossrange  profiles.  As  seen  in  Figure  6.9,  the  amplitudes  of  the 
scatterers'  image  vary  depending  on  the  amount  of  doppler  change  as  seen  by  the 
tracking  radar  during  the  processing  integration  time.  Two  scatterers  located  on  opposite 
sides  of  the  center  of  rotation  will  experience  different  doppler  changes  during  the 
integration  time,  and  consequently  the  amount  of  spreading,  and  hence  their  image 
amplitudes  will  be  different.  The  undistorted  quality  of  the  image  is  shown  in  Figure 
6.10. 
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Figure  6.9  RCS  Returns  (Mesh  Plot)  from  9  Scatterers  on  a  Simple  Target. 
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Figure  6.10  Undistorted  RCS  Returns  (Plan  View)  from  9 
Scatterers  on  a  Simple  Target.  Illustrated  in  Figure  5.1. 
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D. 


IMAGE  OF  COMPLEX  TARGET 


The  image  of  the  complex  target  in  Figure  5.2,  which  simulates  an  aircraft,  is 
plotted  in  Figure  6.11.  The  asymmetry  of  the  image  due  to  the  different  shifting  and 
spreading  of  scatterers  is  evident. 
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Figure  6.1 1  RCS  Returns  (Plan  View)  from  22  Scatterers  on  the  Complex 
Target  (RCS  of  leading  edge  scatterers  higher  than  those  on  trailing  edge). 
Illustrated  in  Figure  5.2. 
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VII.      CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

This  thesis  employed  both  stepped  frequency  and  ISAR  processing  techniques  to 
generate  an  image  from  the  returns  of  targets,  both  simple  and  complex.  The  parameters 
that  affect  the  downrange  and  crossrange  tracking  windows,  range  resolution  and 
ambiguities  were  examined.  Simple  targets  with  few  scatterers  were  used  to  identify  the 
effects  of  radar  parameters  on  image  quality.  This  resulted  in  the  formulation  of 
guidelines  for  an  optimized  window  in  the  design  of  the  imaging  simulation  program. 
Finally,  a  detailed  point  scatterer  model  of  a  complex  aircraft  was  run  through  the 
simulation  and  the  image  plotted.  Scatterers  close  to  one  another  could  easily  be 
discerned  in  the  image. 

The  stepped  frequency  and  ISAR  processing  techniques  are  useful  in  obtaining 
high  resolution  images  of  targets  in  both  downrange  and  crossrange.  Resolutions  of  less 
than  one  meter  are  possible  using  these  two  complementary  techniques.  Together  with 
velocity  compensation  techniques  discussed  in  earlier  thesis  studies  [3],  imaging  of 
targets  comprised  of  multiple  RCS  scatterers  are  possible.  Nevertheless,  it  is  important  to 
choose  the  correct  radar  parameters  in  order  to  avoid  aliasing  and  range  spreading  and 
walking  effects  due  to  the  target's  rotational  motion. 

B.  RECOMMENDATIONS 

Further  investigations  can  be  made  in  the  imaging  techniques.  As  the  ISAR 
processing  is  a  conventional  processing  technique,  it  assumes  that  the  doppler  shifts  are 

51 


constant  during  the  integration  period.  If  the  scatterers  drift  out  of  their  range  cells  or 
their  doppler  frequency  shifts  become  time  varying,  the  image  will  become  blurred. 
Efficient  motion  compensation  techniques  need  to  be  employed  for  each  scatterer,  as 
opposed  to  conventional  motion  compensation  which  applies  compensation  to  the  target 
as  a  whole.  A  recommended  approach  is  to  use  the  Joint  Time  Frequency  Analysis  [6]. 
This  method  replaces  the  conventional  radar  Fourier  transform  with  a  high-resolution 
time-frequency  transform,  and  maps  2-D  range-Doppler  image  frame  into  a  3-D  time- 
range-Doppler  image  cube.  By  sampling  in  time,  a  time  sequence  of  2-D  range-Doppler 
images  with  superior  resolution  can  be  obtained. 
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APPENDIX  A:  STEPPED  FREQUENCY  WAVEFORM  PROGRAM  CODE 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  Program  Code  for  Generation  of  Effective  Bandwidth  of  Step  % 
%  Frequency  Waveform  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear  all 


^Initialize 
format  long 
pw  =  0.1e-6; 
fc=1.0e9; 
Tp  =  4*pw; 
fs  =  1/Tp; 
pul=12; 
del_f=zeros(l,pul); 

for  nn=l:pul 

del_f(nn)=5e6*(nn-l); 
end 


%pulsewidth 

%carrier  frequency 

%PRT 

%sampling  frequency 

%number  of  pulses  per  train/burst 

^initialise  stepped  frequecies 

%stepped  frequency  step  size 


%Compute  f(t)  with  64  sample  points  per  pulse 


datapoint=512; 
del_t=pw/datapoint; 
f=zeros(  1  ,pul  *(Tp/p  w)  *datapoint) ; 
for  pulse=l:pul; 
templ=(pulse-l)*Tp/pw*datapoint; 

for  nn=(templ):((datapoint)+templ) 
n=floor(nn); 

t=n*del_t  +  (del_t/2);  °Ii 

f(n+l)=cos(2*pi*(fc+del_f(pulse))*t); 

end 
end 


%number  of  data  points 
%smallest  time  resolution 
%Initialisse  f(n) 

%define  the  start  of  each  of  pulses 
%(templ=64x4=256  pts,ie,  0,256,512,. 


,32512  or  64x4x127)) 


%  integer  value 

%  take  time  of  each  sub-pulse  to  be  the  mid-point  bet  2  sub-pulses 
%  define  the  amplitude  of  the  pulses 


%Zero  padding  after  last  pulse  for  better  frequency  resolution 

padzero=2*puI*(Tp/pw)*datapoint;%add  128  NULL  pulses  after  the  last  pulse 

last=pul*Tp/pw*datapoint;  %last  point  on  pulse  train 

int=floor(last);  %integer  value 

f((int+l):padzero)=zeros(l,padzero-int);%pad  with  zeros 

A=[0:padzero-1]; 

t(l:padzero)=  del_t*A;  %  set  time  scale 

%Plot  the  pulse  train 

figure(l); 

subplot(2,l,l) 

tlimit  =  pw/datapoint*(padzero-l); 

plot(t,f),xlabeie\fontsize{  8  JTimeCsec)), 

ylabel(^fontsize{  8 }  Amplitude) 

title([\fontsize{8}Pulse  train  padded  with  ',num2str(padzero),  'zeros7]) 

axis([0  tlimit  -1  1]); 

grid  on 
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%Compute  spectrum  of  pulse  train 

subplot(2,l,2) 

FT=((fft(f))); 

temp=length(abs(FT)); 

plot(abs(FT)),xlabel(\fontsize{8}Bin  number),ylabel(\fontsize{  8}  Amplitude) 

title(\fontsize{8}FFT) 

grid  on 

%Plot  spectrum  using  FFTSHIFT 

figure(2) 

subplot(2,l,l) 

FFTSHT=fftshift(FT); 

plotCabsCFFTSHT^^labeKXfontsizefSJBinnumber^ylabeKNfontsizefS}  Amplitude') 
titIe(\fontsize{8}FFTSHIFT) 
grid  on 

%Plot  spectrum  in  frequency  domain 
subplot(2,l,2) 

deltafreq  =  l/(pw/datapoint*padzero);  %  frequency  stepsize 

f  =  ((l:padzero)-(padzero/2))*deltafreq; 

plot(f,abs(FFTSHT)); 

x!abel(\fontsize{  8  {Frequency  (Hz)),ylabel(\fontsize{  8 }  Amplitude) 

grid  on 

Fmax=max(abs(FFTSHT)) ; 

figure(3) 

deltafreq  =  l/(pw/datapoint*padzero); 

f  =  ((l:padzero)-(padzero/2))*deltafreq; 

ff=abs(f); 

pIot(f,abs(FFTSHT)/Fmax); 

xlabelCyfontsize{8}Freq  in  Hz ),ylabel(\fontsize{  8}  Amplitude) 

grid  on 

figure(3),  axis([0.9e9  1.15e9  0  1]) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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APPENDIX  B  :  STEPPED  FREQUENCY  AND  ISAR  PROCESSING  PROGRAM 

CODES 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  Program  Code  for  Stepped  Frequency  and  ISAR  Simulation  in  Imaging  Radar  % 

%%%%%%%%%%%%7c%%%%%%%%%%%%%%%%%%%%%%%%%% 

^Multiple  scatterers  rotating  about  a  reference  point 

%Process  a  pulse  train  for  doppler  and  range  information 
clear  all 
format  long 

%%%%%%%%%%% 

^INITIALIZATIONS 

%  %  %  %  %  %  %  %  %  %  % 

%Radar  parameters 

N=5 12;  %  number  of  pulses  in  one  burst 

fp=10e3;  %  PRF 

Tp=l/fp;  %  interpulse  period 

T=N*Tp;  %  total  observation  time 

df=l/T;  %  frequency  increment 

fc=40e9;  %  carrier  frequency 

c=3e8;  %  speed  of  light 

wave=c/fc;  %  wavelength 

step=le6;  %  stepped  frequency  step  size 

taucom=l/(N*step) ;  %  compressed  pulse  width 

dr=c*taucom/2;  %  compressed  down-range  resolution 

fpeff=fp/N;  %  effective  PRF  using  stepped  frequency  (used  in  cross  range) 

%Define  stepped  frequency 
for  nn=l:N 

del_f(nn)=step*(nn-l); 
end 


%%%%%%%%%% 
%Target  generation^ 
%%%%%%%%%% 


%Define  Radar  position 

radarpos=[0  0  0] ;  %define  radar  position 

%Define  initial  scatterers' positions  on  target 
x=[4530  4525  4527  4527  4530  4530  4540  4540  4540]; 
y=[0      0    -3.5  3.5    -7    7     0    -1.5   1.5]; 
z=[0      0     0     0     0    0     0     0     0]; 

center=[4530  0  0];  %center  of  rotation 

r  =  sqrt((x-radarpos(l)).A2+(y-radarpos(2)).A2+(z-radarpos(3)).A2); 
%Radial  distance  of  scatterers  from  radar 

%Compute  no.  of  down-range  bins 
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Ru=c/(2*step);  %unambiguous  down-range  window  =  150m  (0.29mx512) 

Rumax=c/(2*fp);  %unambiguous  maximum  detection  range  (15000m) 

numdrbin=  Rumax/Ru;       %no.  of  down-range  bins***(  15000/150=100) 

%Compute  number  of  pulse  burst 

temp=mod(r,Ru)/Ru; 

mbin=temp*N; 

rngno      =  floor(r./Ru);       %  range  gate  no. 

ntarg=length(mbin);  %  no  of  scatterers 

M=floor(Rumax/Ru) ; 

%  Assign  velocity  to  each  scatterer 

w=0.002;  %  rotation  speed  of  target  in  rad/s  in  x-y  plane*** 

%positive  =  clockwise 

%Determine  doppler  velocity  for  each  scatterer  at  each  time  frame,  ie,  every  N  pulses  for  M  bursts. 

ANS=  angle2(w,center,x,y,z,radarpos,ntarg,M,N,Tp,wave);    %doppler  velocity 

fd=ANS(:,l:M); 

vr=ANS(:,M+l:2*M); 

range_time=ANS(:,2*M+l:3*M);  %range  of  scatterer  at  different  time  throughout  Ti 

disp([No.  of  scatterers  defined  :  \num2str(ntarg)]); 
dispCfMax.  unambiguous  detection  range  :  \num2str(Rumax)]); 
disp([ Unambiguous  range  window  :  ',num2str(Ru)]); 
disp([The  calculated  coarse  downrange  bin  #  :',num2str(rngno)]); 
disp([  Calculated  fine  downrange  bin  #  :',num2str(mbin)]); 

%Generate  Is  and  Qs 

I=zeros(N,M); 

Q=zeros(N,M); 

dt=2*dr/c;  %time  separation  between  samples 

A=ones(l,ntarg); 
B=5*ones(l,ll); 
A(1:11)=B;  %define  dominant  scatterers  on  leading  edge  of  aircraft 

for  tgt=l:ntarg 

for  time=l:M*N  ; 

tm(tgt,time)=(2/3e8)*(range_time(tgt,ceil(time/N))  );  %2Rn/c 
end 

for  m=l:M 
forn=l:N 
I(n,m,tgt)=A(tgt)*cos(2*pi*(fc+del_f(n)+fd(tgt,m)).*tm(tgt,((m-l)*N+n))); 
Q(n,m,tgt)=A(tgt)*sin(2*pi*(fc+del_f(n)+fd(tgt,m)).*tm(tgt,((m-l)*N+n))); 
end 
end 
end 

II=zeros(N,M); 

QQ=zeros(N,M); 

for  tg=l:ntarg 

II=II+I(:,:,tg); 

QQ=QQ+Q(:,:,tg); 

end 
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%Aadd  noise  if  required 

No=randn(N,M)*0.05*0; 

I=H+No;  Q=QQ+No; 

%save  c:\files\echoes  I  Q  -ASCII 

%Plot  the  position  of  scatterers  wrt  radar 

figure(l)  %3-D  plot 

subplot(2,l,l) 

plot3(x,  y,  z,*1") 

xlabelfx  co-ordinates') 

ylabel('y  co-ordinates') 

zlabel(  z  co-ordinares") 

grid  on;  view(3); 

figure(  1 ),  hold  on  %Plan  view 

subplot(2,l,2) 

plot3(x,  y,  z,"<") 

xlabel('x  co-ordinates') 

ylabel('y  co-ordinates') 

zlabel('z  co-ordinares) 

grid  on;  view(45,90); 

axis([min(x)-5,  max(x)+5,  min(y)-5,  max(y)+5,  -1,  1]) 

%%%%%%%%%%%%%%%%%%%%%% 

%processing  of  target  returns% 

%%%%%%%%%%%%%%%%%%%%%% 

%Processing  Target  Returns  using  IFFT  to  obtainb  HRR  profile 

it=0;  %initialise  no.  of  scatterers  detected 

bincourse=0;  binfine=0;  dopbin=0;  %initialise  range  and  doppler  bin  of  detected  scatterer 

Zl=zeros(M,N); 

%FFT  of  the  time  samples  will  give  the  HRR  profile 
form=l:M  %M  range  bins 

Al=I(:,m)-j*Q(:,m); 

Bl=fft(Al);  Asl=ifft(Al); 

Zl(m,:)=Asl(:).'; 
end 

figure(2)  %plot  HRR 

rngfaxis=([0:N-l]*dr); 

fort=l:M 

plot(rngfaxis,abs(Zl(t,:)))  %/N) 

hold  on 
end 

xlabel(Downrange  in  meters') 
ylabel(  Amplitude') 
grid  on 

%Processing  Target  Returns  using  FFT  to  obtain  crossrange  profile 

Z2=zeros(N,M); 
forp=l:N 
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B2=fft(Zl(:,p));B3=fftshift(B2); 
Z2(p,:)=B3(:).'; 

end 

figure(3)  %plot  doppler  profile 

dfstep=fp/(N*M); 

drc=(3e8*dfstep)/(2*w*fc); 

craxis=(([0:M-l]-M/2)*drc); 

forn=l:N 

plot(craxis,abs(Z2(n,:))/M) 

hold  on 
end 

xlabel('Crossrange  in  meters') 
ylabelf Amplitude) 

grid  on 


figure(4) 

x=linspace(  -((M-l)-M/2)*drc 

y=linspace(0,  (N-l)*dr,  N); 

mesh(x,y,abs(Z2)) 

xlabeK'Crossrange  in  meters') 

ylabel(Downrange  in  meters) 

zlabelCAmplitude7) 


,((M-l)-M/2)*drc,  M); 


figure(5) 

x=linspace(  -((M-l)-M/2)*drc  ,((M-l)-M/2)*drc, 

y=linspace(0,  (N-l)*dr,  N); 

mesh(x,y,abs(Z2)) 

xlimit=((M- 1  )-M/2)  *dfstep; 

ylimit=(N-l)*dr; 

zlimit=0; 

M); 

forf=l:M 

zlimit(0=max(abs(Z2(:,f))); 
end 

zzlimit=max 
n=21; 

(zlimit); 

mymap=[  1.0    1.0    1.0; 
1.0   1.0  0.9; 

1.0 

1.0  0.8; 

1.0 

1.0   0.7 

1.0 

1.0    0.6 

1.0 

1.0    0.5 

1.0 

1.0    0.4 

1.0 

1.0    0.3 

1.0 

1.0    0.2 

1.0 

1.0  0.1; 

1.0 

1.0  0.0; 

1.0 

0.9  0.0; 

1.0 

0.8  0.0; 

1.0 

0.7    0.0 

1.0 

0.6    0.0 

1.0 

0.5    0.0 

1.0 

0.4    0.0 

1.0 

0.3    0.0 
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1.0  0.2  0.0; 
1.0  0.1  0.0; 
1.0    0.0    0.0;]; 

colormap(mymap) 

imagesc(x,y,abs(Z2)); 

colorbar('hori') 

figure(5) 

xlabel('Crossrange  in  meters') 

ylabel(Downrange  in  meters') 

zlabel( 'Amplitude') 

view(0,90) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%Function  to  Compute  Radial  Velocity  of  Individual  Scatterers  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  ANS=angle2(w,center,x,y,z,radarpos,ntarg,M,N,Tp,wave) 

format  long 

radius=sqrt((center(2)-y).A2+(center(l)-x).A2);     %radius  of  scatterer  from  center  of  rotating  circle,  vector 

range  =  sqrt((x-radarpos(l)).A2+(y-radarpos(2)).A2+(z-radarpos(3)).A2%center-radar  distance,  vector 

range_ref=sqrt((center(l)-radarpos(l)).A2+(center(2)-radarpos(2)).A2+(center(3)-radarpos(3)).A2); 

v=w*radius; 

%Determine  sign  of  radial  velocity  from  gradient  formula 

delta_x=x-radarpos(l); 

delta_y=y-radarpos(2); 

grad_ref=(center(2)-radarpos(2))/(center(l)-radarpos(l));  %gradient  of  center-radar  vector 

grad=(delta_y)  ./(delta_x) ; 

%Determine  x  and 
c=range; 
b=radius; 
a=range_ref; 

%  Vector  OR  :  radar-center  vector 

OR=[center(  1  )-radarpos(  1 ) ;  center(2)-radarpos(2)] ; 

%Vector  RA  :  center  to  scatterer 

for  tg=l:ntarg 

RAx(tg)=[x(tg)-center(l)]; 

RAy(tg)=[;y(tg)-center(2)]; 

end 

RA=[RAx;RAy]; 


%Compute  length  of  OR  and  RA 

ORmag=range_ref; 

RAmag=radius; 
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%Compute  angle  between  OR  and  RA  ie  radar-center  and  center-scatterer 
%  Quadrant  1 

for  tg=l:ntarg 
if  RAmag(tg)~=0 
thetal(tg)=acos(  dot(OR,RA(:,tg))  /(ORmag*RAmag(tg))  );  %dot  product  between  vectors  OR  &  RA 
%verified 
else 

thetal(tg)=pi; 
end 
end 

%Compute  theta  for  M  pulses  &  x  and  y  co-ordinates  of  scatterers  over  time 
for  tg=l:ntarg 

%check  which  quadrant  scatterer  is  in 
if  ( (  grad(tg)>grad_ref )  |  ((  grad(tg)==grad_ref )  &  (range(tg)>range_ref))  )  %quadrant  1 

for  m=l:M 

theta(tg,m)=theta  1  (tg)+w*N*Tp*(m- 1 ); 
x(tg,m)=radius(tg)*cos(theta(tg,m)); 
y(tg,m)=radius(tg)*sin(theta(tg,m)); 
if  m==l 

range_time(tg,m)=range(tg);     %initial  range  of  scatterer  before  rotation 
else 

range_time(tg,m)=sqrt(aA2+b(tg)A2-2*b(tg)*a*cos(pi-theta(tg,m)) ); 
%range  of  scatterer  over  time 
end 
%discard  last  value 
xcoord(tg,m)=center(l)+x(tg,m);  %x-co-ordinates  of  A  %diff  sign  for  quad  2 

ycoord(tg,m)=center(2)+y(tg,m);  %y-co-ordinates  of  A 

OA=[xcoord(tg,m)-radarpos(l),  ycoord(tg,m)-radarpos(2)];  %  vector  OA 

RA=[xcoord(tg,m)-center(l),ycoord(tg,m)-center(2)];  %magnitude  of  RA=radius(tg) 

OAmag(tg,m)=  sqrt((xcoord(tg,m)-radarpos(l))A2  +  (ycoord(tg,m)-radarpos(2))A2  );  %magnitude  of 
OA 

%angle  bet  dot  product  of  vectors  OA  &  RA 
beta(tg,m)=acos(  dot(OA,RA  )/(OAmag(tg,m)*RAmag(tg))); 
del_theta(tg,m)=abs(  pi/2-beta(tg,m) ); 
vr(tg,m)=-v(tg)*cos(del_theta(tg,m)); 
fd(tg,m)=-2*vr(tg,m)/wave; 

end 
elseif  ( (  grad(tg)<grad_ref )  |  ((  grad(tg)==grad_ref )  &  (range(tg)<range_ref))  )  %quadrant  2 

for  m=l:M 

theta(tg,m)=thetal(tg)-w*N*Tp*(m-l); 
x(tg,m)=radius(tg)*cos(theta(tg,m)); 
y(tg,m)=radius(tg)*sin(theta(tg,m)); 

ifm==l 

range_time(tg,m)=range(tg);     %initial  range  of  scatterer  before  rotation 

else 

range_time(tg,m)=sqrt(aA2+b(tg)A2-2*b(tg)*a*cos(pi-theta(tg,m))  );  %range  of  scatterer  over  time 

end 

%Discard  last  value 

xcoord(tg,m)=center(l)+x(tg,m);  %x-co-ordinates  of  A  %diff  sign  for  quad  2 

ycoord(tg,m)=center(2)-y(tg,m);       %y-co-ordinates  of  A  %diff  in  sign 
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OA=[xcoord(tg,m)-radarpos(l),  ycoord(tg,m)-radarpos(2)];      %  vector  OA 
RA=[xcoord(tg,m)-center(l),ycoord(tg,m)-center(2)];  %magnitude  of  RA=radius(tg) 

OAmag(tg.m)=  sqrt((xcoord(tg,m)-radarpos(l))A2  +  (ycoord(tg,m)-radarpos(2))A2  ); 
%magnitude  of  OA 

beta(tg,m)=acos(  dot(OA,RA)/(OAmag(tg,m)*RAmag(tg))     ); 
%angle  bet  dot  product  of  vectors  OA  &  RA 
del_theta(tg,m)=abs(  pi/2-beta(tg,m)); 
vr(tg,rn)=v(tg)*cos(del_theta(tg,m)); 
fd(tg,m)=-2  *  vr(tg,m)/wave; 
end 

elseif(  (grad(tg)==grad_reO  &  range(tg)==range_ref )  %scatterer  is  at  center 
for  m=l:M 
theta(tg,m)=pi; 

x(tg,m)=radius(tg)*cos(theta(tg,m)); 
y(tg,m)=radius(tg)*sin(theta(tg,m)); 
ifm==l 

range_time(tg,m)=range(tg);     %initial  range  of  scatterer  before  rotation 
else 

range_time(tg,m)=sqrt(aA2+b(tg)A2-2*b(tg)*a*cos(pi-theta(tg,m)) ); 
%range  of  scatterer  over  time 
end 

end 


end 
end 


ANS=[fd,vr,range_time(:,l:M)]; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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